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Abstract 

The Chern numbers which correspond to quantized Hall conductance <x A3 , were calculated for single- and bi-layer honeycomb lat- 
tices. The quantization of cr xy occurs in entire energy range. Several large jumps of Chern numbers appear at van-Hove singularities 
of energy bands without magnetic fields. The plateaux of cr xy are discussed from semi-classical quantization. 
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1. Introduction 

The quantum Hall effect (QHE) is one of the most peculiar 
phenomena in two-dimensional electron system. The quantized 
Hall conductance <x„, has been widely accepted as a topological 

quantity HUSHIlSHIl]. When the chemical potential is 
located within an energy gap, the cr A ,, is given as 



CTjv = — -c F , 
h 



(1) 



where c F is a Chern number, which is an topological integer de- 
fined for occupied states. The c F has been computed for simple 
tight-binding models and exotic nature of electron states under 
uniform magnetic field has been revealed. 

Recently, the discovery of anomalous QHE in graphene |0, 
[loll shed a new light to the quantization of Hall conductance. 
The <x n . of graphene shows quantized plateaux at cr xy = {2N + 
l)e 2 /h per spin, where N is an integer. This means that only 
plateaux with odd number appear and the plateau corresponding 
to a xy = does not occur. This anomalous QHE originates 
from the massless Dirac particles realized at the K and K' points 
in the Brillouin zone of graphene. 

In our previous study 111 ill , we calculated Chern numbers 
for realistic multi-band tight-binding models, which describe 
electrons on graphene. We found that anomalous quantization 
of Hall conductance persists up to van-Hove singularities, as 
shown in the previous calculations using simplified model [12]. 
In addition, the envelope of quantized Hall conductance can 
be well explained with semi-classical expressions. We further 
found that the Onsager's quantization rule IU3I1 . 



Si(e) 

Q R 7 



= (n + y)(j), 



(2) 
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can predict the positions of plateaux when the band dispersion 
is well described by single band. Here, <p is magnetic flux pene- 
trating in a unit cell and Qbz is the area of irreducible Brillouin 
zone. For Dirac particle regions, above equation also repro- 
duces the positions of Landau levels if we choose Maslov index 
y — 0. As this expression does not rely on the linear dispersion 
of Dirac particles, it can be applicable in wide energy region. 
Therefore, this expression may be utilized for analysis of QHE 
and other phenomena of Landau quantization, e. g., de Haas- 
van Alphen effects. 

In this paper, we extend our study to modified version of 
tight-binding models related to graphene layers. We study how 
the sequence of Chern numbers may be altered and whether 
semi-classical rules can be applicable to the extended models. 

2. Calculation Methods 

We describe electrons on honeycomb lattice by tight- 
binding models. The spin degeneracy is ignored for sim- 
plicity. The uniform magnetic field is introduced as phase 
factors 6jj of hoppings so that they satisfy Zciosedioop Qij — 
27r(flux quanta in the loop). 

As the original translational symmetry is broken under uni- 
form magnetic field, We must consider magnetic translational 
symmetry, which is characterized by the magnetic flux <p = 
BCl/ipu where Q is the area of a unit cell and ipo = hc/e is the 
quantized magnetic flux. When is a rational number p/q,we 
recover enlarged magnetic unit cell whose area is q times larger 
than that of the original cell. Then, eigenstates are labeled by a 
wave-number k. Thus, an energy band splits into q sub-bands 
by the uniform field. 

Although the Chern number Cf is given by a sum of the 
Chern numbers of all filled bands [1], evaluation of them sepa- 
rately causes numerical difficulty since one needs to treat many 
Chern numbers for the weak field limit, where most of them 
cancel with each other. Physical Chern number of the filled 
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Figure 1: (a) A honeycomb lattice with staggered on-site energies. Two in- 
equivalent sites A and B are distinguished by white and black circles. The 
dashed lines indicate hoppings between next-nearest sites, (b) A structure 
model of bilayer graphene. Carbon atoms on upper and lower layer are in- 
dicated by small and large circles, respectively. 



many body state is stable if the energy gap above the Fermi en- 
ergy is stable. It is independent of a possible level crossing at 
far below the Fermi energy. Then the non Abelian formulation 
of the Chern number by collectin g all of the filled bands into a 
single multiplet is quite useful Il4ll5lll6ll . Then the Chern 
number Cf of the multiplet ^ is calculated by discretizing the 
Brillouin zone into mesh {k(} as II 1 611 



C F (jl) 

F(k) 
u^k) 



2 -V 



Arg u x (k)u y {k + Ak x )[u x (k + Ak y )u y (k)]* 
det[ v F(£)] tv F(A: + Akp) 

(\Mk)),--- ,liM*)» 



where Ak v is a discretized momentum along v direction, 
\ift„,(k)) is a Bloch state of the band index m and m = M is 
the highest occupied energy band below the chemical poten- 
tial p. This expression can be understood as a two-dimensional 
analogue of the KSV formula II 1711 . 

This formulation was successfully applied to quantum Hall 
effect on a disorder svstem lfl8ll . graphene lfl2l . When the field 
<p = p/qis sufficiently small, p sub-bands are grouped together 
in general. These grouped bands correspond to a Landau band 
of nearly free electrons. Hereafter, we treat the field <p — 1/q, 
which enables us to consider the weak-field limit easily. 



3. Honeycomb lattice with staggered on-site energy 

We treat a honeycomb lattice with two inequivalent sites A 
and B as shown in Fig.[TJa). Single ^--orbital is put on each site. 
The hoppings between nearest and next-nearest neighbor sites 
are included with parameters t\ and t<x. The on-site energies are 
set to ±A for A and B sites. Hereafter, we set t\ = 1, i.e., t\ is 
chosen as an energy scale. 

When ?2 = A = 0, this model describes simplified 7r-bands on 
graphene and the band dispersions which resemble zero-mass 
Dirac cones exist at two symmetrically equivalent points, K 
(1/3 1/3) and K' (2/3 2/3), in the irreducible Brillouin zone. 
The QHE of this model has been extensively studied theoret- 
ically 111 211 . Here, we briefly explain the correspondence be- 



tween quantum mechanical treatment by Chern numbers and 
semi-classical quantization rule ©. 

The calculated Chern numbers for <p = 1/200 are plotted as 
a function of chemical potential in lower panel of Fig. [2] The 
upper panel shows the total density of states (DOS) when mag- 
netic field is absent. From these figures, we can easily recognize 
that large discontinuous jumps of Chern numbers correspond to 
van-Hove singularities at s = ±1. Within semi-classical theory, 
this discontinuity can be explained as a topological transitions 
of Fermi surface from electron (hole) pockets to hole (electron) 
pocketslflTll. 

Fig. [3] presents enlarged views of Chern numbers near the 
van-Hove singularity at s — 1. Anomalous QHE survives up 
to van-Hove singularities and the sequence of Chern numbers 
is odd number only, while above singularities Chern numbers 
change with steps of 1, recovering ordinary QHE 11211 . For 
both regions, semi-classical quantization rule can explain 
the positions of Landau levels if we choose y — 1 /2 at ordinary 
QHE region and y — at anomalous QHE region. This result 
is non-trivial near the van-Hove singularity as the band disper- 
sion around s = 1 is significantly deformed from the continuum 
expressions of massless Dirac s ~ \k\ or normal effective mass 
theory j-?. 

We extend our calculations by incorporating next nearest 
neighbor ?2 and staggered on-site energy ±A. When A = 0, the 
hoppings ?2 alone do not destroy the massless Dirac cones. In 
this case, the calculated Chern number sequence remains with 
cf = ■ ■ • , -5, -3, -1, 1, 3, 5, • ■ ■ as shown in Fig. 0a). The 
positions of Landau levels are still explained with the semi- 
classical quantization with y — 0. 

With finite A, the Dirac cones at K and K' acquire band 
gaps of order 2A, which means that Dirac particles have fi- 
nite "mass". When the chemical potential is located above the 
gap, Fermi surface is still constructed by two electron pock- 
ets around K and K' . In addition, the electronic states around 
K and K' are related by lattice symmetry. If we interpret this 
band structure naively, we may expect that two ordinary elec- 
tron pockets generate Landau levels at same energies and that 
Chern numbers between adjacent energy gaps differ by 2. How- 
ever, we found that such interpretation does not hold for current 
model. As plotted in Fig.|4jb), additional plateaux appear in cr xy 
and the Chern numbers increase with steps of 1. In addition, 
the positions of Landau levels do not agree with semi-classical 
rule© with y = nor 1/2. These results mean that electron 
pockets near K and K' generate Landau levels at different en- 
ergy positions even though they are equivalent when magnetic 
field is absent. Indeed, the positions of Landau levels may be 
explained if we choose different y for Fermi surface segments 
around K and K' . 



4. Bilayer graphene 

The graphene bilayer consists of two honeycomb lattices as 
shown in Fig. [Tib). Each layer is composed of two inequivalent 
sites which we denote A \ and B \ on lower layer, and A 2 and B2 
on upper one. The B2 sites on upper layer are located above 
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Figure 2: The upper panel shows the total density of states (DOS) without 
magnetic field. The lower panel presents calculated Chem number sequence as 
a function of chemical potential. 
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Figure 3: The upper panel shows the total density of states (DOS) without 
magnetic field. The lower panel presents calculated Chern number sequence as 
a function of chemical potential. The vertical dashed lines indicate the positions 
where semi-classical quantization rule is satisfied with (a) y = and (b) y = 
1/2. 
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Figure 4: The total DOS and Chern numbers for a honeycomb lattice model. 
The magnetic field is (p = 1/101. The dashed lines indicate the positions where 
semi-classical quantization rule is satisfied with y = 0. (a) A = 0, t% = 0.3. (b) 
A = 0.5, t 2 = 0.3. 

the A\ sites on lower layer. The largest hoppings yo between n- 
orbitals are between nearest-neighbor atoms on the same layer. 
The inter-layer hoppings are included between A\ and Bo sites 
with parameter j\ . Other hoppings are ignored for simplicity. 
These parameters are taken from previous estimations II 1911 as 
yo = 2.6eV and j\ /yo ~ 0.14. 

The OHE of bilayer graphene was observed experimen- 
tally 12011 and its electronic structure has been studied 
ex tensivelyj23|2l|2l. Within the present model, the electron- 
hole symmetry is maintained and band dispersions are sym- 
metric around s = 0. The inter-layer hoppings modify the 
Dirac cones which exist in single-graphene and results in two 
branches of energy bands for s > 0. Both branches behave 
as s — 6k 2 near the K and K' points. In addition, the lower 
branches remain gapless. 

The uniform magnetic field is applied along perpendicular 
to layers. Fig. [5] shows global structure of calculated Chern 
numbers as a function of chemical potential. The total DOS 
without magnetic field is also plotted with same energy scale. 
Because two branches of energy bands exist, pair of van-Hove 
singularities appear around s — ±1. Accordingly, the large dis- 
continuous jumps of Chern numbers around s — +1 split into 
two jumps. At s — 0, the Chern numbers jump from -2 to 2, 
suggesting the existence of 4 Landau levels at e = 0. The po- 
sitions of other Landau levels shown in Fig. |6ja) seem to be 
located without concrete rule. This is because both upper and 
lower branches of energy bands contribute to the quantized Lan- 
dau levels. Nonetheless, the Chern numbers always jump by 2, 
which results in the sequence of cp = ■ ■ ■ , -4, -2, 2, 4, • • • . 

We studied Landau levels from lower branch of energy bands 
by lifting the upper branch with large j\. Fig. |6jb) shows the 
sequence of Chern numbers with yi /yo = 0.6. With this param- 
eter, only lower branch with gapless dispersion contributes to 
quantized Landau levels near s = 0. As a result, the appearance 
of Landau levels becomes systematic but can not be explained 
by standard semi-classical rule irrespective to the choice of y. 
In the continuum limit, analytical expression of Landau levels 
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position of Landau level if we choose y appropriately. When 
the on-site energies are introduced, additional plateaux will ap- 
pear. This indicate that valley degeneracy between K and K' is 
broken when the magnetic field is applied. 

For bilayer graphene, the Chern number sequence of Cf = 
• • • , -4, -2, 2, 4, 6, ■ • • is obtained, which agrees with previous 
theoretical results. The semi-classical interpretation is difficult 
when multiple bands contribute to the same energy regions. 
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Figure 5: The total DOS (upper panel) and Chern number sequences for bilayer 
graphene are shown. 
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Figure 6: The total DOS and Chern numbers for bilayer honeycomb lattice. 
Parameters are chosen as (a) yi/yo = 0.14 and (b) yi/yo = 0.6. The vertical 
lines indicate the positions of Landau levels calculated from modified semi- 
classical rule 



has been obtained as s ~ yln(n - 1) [21]. Using this expres- 
sion as a guide, we have tried several modifications to the semi- 
classical rule and found that the equation 
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agrees with calculated Landau levels within the energy scale 
presented in Fig. 0b). However, we have not found any the- 
oretical support for above equation. Further studies would be 
necessary to clarify the quantization rule. 



5. Conclusion 

We studied quantum Hall effect on single and bilayer honey- 
comb lattices. For single-layer honeycomb lattice without stag- 
gered on-site energies, the semi-classical theory explains the 
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